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Abstract. We show by detailed model computation how much the Pre-Main Sequence (PMS) Lithium depletion depends on 
the treatment of over-adiabaticity, by taking advantage of the results of new models by Montalban et al., which apply different 
treatments of convection to non-grey PMS models. In order to reproduce both the PMS lithium depletion (inferred from the 
lithium depletion patterns in young open clusters), and the location of PMS tracks in the HR diagram (inferred from the study 
of young PMS stars), convection both in the atmosphere and in a good fraction of the stellar envelope must be highly inefficient: 
e.g., in the Mixing Length Theory approximation, it must have a very low a = l/H p . Unfortunately, the radii of these models are 
at variance with the solar radius, possibly indicating that there is some additional physical input, generally not taken into account 
in the stellar models, which affects the efficiency of convection in PMS stars, but probably not in the main sequence stars nor 
in evolved red giants. We stress the importance of determining precisely masses and lithium abundance in PMS binaries such 
as the important spectroscopic and eclipsing binary RXJ 0529.4+0041. 
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1. Introduction 

The location in the HR diagram of the Pre-Main Sequence 
(PMS) evolutionary tracks is very sensitive to physical inputs 
such as low-temperature opacity, equation of state, rotation, at- 
mosphere model, and the treatment of convection. During the 
last years, much work has been done to improve the knowledge 
of low-temperature opacities and to include them in the model- 
ing of stellar atmospheres. Among the many theoretical evolu- 
tionary tracks today available (new convection treatment with 
grey boundary conditions; classic convection with non-grey at- 
mosphere models...) it is difficult to disentangle the effect of the 
different physical inputs on the results. In paper I (Montalban 
et al. 2003) we made a detailed exploration of models, with 
the aim to extricate the different roles of non-grey atmospheres 
and of convection in the track computations. This task has been 
approached by comparing a large number of models, having 
different assumptions for the treatment of convection both in 
the atmosphere and in the interior, and using different grids of 
atmospheric computations. In particular, we have shown that 
a stellar model is fully described only when we specify not 
only the atmospheric model used as boundary condition to the 
interior, but also the convection parameters used for the atmo- 
spheric grid, and the value of optical depth at which the bound- 
ary conditions are taken, that is the photospheric 'matching 
point' T p h. In fact, until now, the 'expensive' model atmosphere 
computations have been generally performed only for one spe- 
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cific convection model, e.g. a given ratio of mixing length to 
pressure scale heigth, in the Mixing Length Theory (MLT), 
a = a a tm- The value of a can be changed only in the com- 
putation of the interior (a = a m ) e.g. to fit the solar radius in 
the solar model. However, if a large value of r p h is chosen as 
matching point between the atmosphere and the interior, the 
most superadiabatic part of the convection zone is all included 
in the atmosphere. Consequently, the changing of a m does not 
affect the model in the same way as a full change the convec- 
tion parameter in the whole model, including the atmosphere 
(see paper I). 



Although this may seem only of academic interest, it is of 
subtle importance, if we wish to progress in one of the most in- 
triguing problems of the stellar structure, namely the problem 
of PMS lithium depletion. In fact, the well known 'problem of 
lithium in the Sun', which in the years 1965-1990 was mainly 
that solar models could not burn any substantial fraction of their 
initial lithium during the PMS, was generally taken as a good 
proof that additional mechanisms for depletion were required, 
acting during the long solar MS lifetime. This interpretation is 
today taken as the most plausible one, confirmed by the varia- 
tion of the Lithium Depletion Patterns (LDP) in open clusters 
(see e.g. Chaboyer 1998). In fact the lithium vs. T e s relation for 
the MS stars of young open clusters indicates a lithium deple- 
tion by at most a factor two for the solar mass in young clusters, 
while it is compatible with the solar depletion (a factor ~140 
with respect to the solar system abundance) in some stars of the 
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Table 1. Lithium remnant after the PMS depletion in 1M ( 
models, starting with log N(Li); n = 3.31 



atmosphere 


Q'atm 


interior 




Tpb 


log N(Li) 


ATLAS9 MLT 


0.5 


MLT 


1.85 


3 


1.864 


ATLAS9 MLT 


0.5 


MLT 


2.3 


100 


2.130 


ATLAS9 MLT 


0.5 


MLT 


6.3 


10 


2.564 


ATLAS9 FST 


0.1 


FST 


0.2 


10 


1.537 


AH97 MLT 


1.0 


MLT 


1.9 


3 


1.960 


AH97 MLT 


1.0 


MLT 


1.9 


100 


2.62 



cluster M67, close to the solar age. For recent reviews see e.g. 
Jeffries (2000) and Pasquini (2000). 

However, a different problem emerges from the most recent 
computation of solar models: they deplete too much lithium 
during the PMS evolution (D' Antona and Mazzitelli 1994, 
1997, Schlattl and Weiss 1999, Piau and Turck-Chieze 2002) 
and are incompatible with the open clusters observations. 1 This 
problem is most severe in models using very efficient convec- 
tion models, in fact it is more relevant for the D'Antona and 
Mazzitelli (1994 and 1997) models adopting the Full Spectrum 
Turbulence (FST) convection by Canuto and Mazzitelli (1991) 
and/or by Canuto et al. (1996). MLT models of the most recent 
generation, adopting updated equations of state and opacities 
also deplete too much lithium: notice however that the problem 
is found mainly in tracks whose convection is adjusted to pro- 
vide the solar radius at the solar age! If one does not require the 
solar fit, it is easy to decrease the convection efficiency and to 
obtain a small PMS lithium depletion. 

A similar problem is found when comparing the location 
of PMS theoretical tracks with the observed few data of PMS 
stars for which an independent determination of mass is avail- 
able, either because they belong to binaries (e.g. Covino et al. 
2001, Steffen et al. 2001) or by the measure of the dynami- 
cal properties of their protoplanetary disks (Simon et al. 2000): 
the tracks most consistent with the observations are those with 
cooler atmospheres (higher mass for a given spectral type) and 
thus those which, generally, provide a radius larger than R Q 
for the solar model. In particular, we can quote the analysis 
by Simon et al. (2000): their results are consistent, -within the 
statistical errors- with the set of tracks whose evolution for the 
solar mass does not pass through the solar location. In addi- 
tion, their observations are not consistent with the FST tracks. 
In fact, for any value of its fine tuning parameters, the FST con- 
vection always provides solar tracks which pass within ~ 100K 
of the solar location. 

These two problems, HR diagram location of the tracks 
during the PMS evolution, and lithium depletion, are neces- 
sarily correlated each other, because, the smaller the T e ff of the 
Hayashi track, the smaller will be the temperature at the base 
of the convective envelope during the possible phase of lithium 
burning (see, e.g. D'Antona et al. 2000). Having computed sev- 
eral sets of stellar models, which are fully described in paper 
I, we decided to look at the lithium depletion with different as- 

1 We mention, as a caveat, that these models are not strictly incom- 
patible with the present solar lithium abundance, if we accept that all 
or most of the solar lithium was depleted at the PMS stage. 



sumptions. In Sect. 2 we show and discuss the results. In Sect. 
3 we show that the models having low lithium depletion also 
reproduce better the location of PMS stars in the HR diagram, 
but they do not reproduce the solar radius. Thus the problem is 
now better defined: the efficiency of convection in PMS must be 
highly reduced with respect to the efficiency in MS. In particu- 
lar, those tracks whose convection formulation can not provide 
a solar model with radius far different from the solar radius 
-noticeably the FST models- can not reproduce the location 
in PMS and its lithium depletion, unless the efficiency of their 
convective transport is altered by another physical input acting 
during the PMS, but not during the MS or post-MS giant phase. 



2. Standard models: fitting the solar radius, or not 

Our models are computed following standard Henyey type inte- 
gration, by means of the ATON2.0 code (Ventura et al. 1998a, 
Montalban et al. 2001) updated in paper I to include as sur- 
face boundary conditions (BCs) the new grids of ATLAS9 at- 
mospheres from Heiter et al. (2002a). We consider here the 
models having solar metal abundance in mass fraction Z=0.02 
and helium abundance Y=0.28. The initial lithium abundance 
was taken as log N(Li)=3.31±0.04, the solar system abundance 
given by Anders and Grevesse 1989. The PMS lithium burning 
starts when the temperature reaches ~ 2x 10 6 K during the PMS 
contraction. The important inputs are the temperature and den- 
sity at the bottom of the convective region. Actually lithium 
burning generally begins approximately at the bottom of the 
Hayashi track, when a radiative core develops and grows dur- 
ing the evolution towards the main sequence. Figure ^ shows 
the the central density, the density at the bottom of the convec- 
tive region, and the surface lithium abundance along two tracks 
of 1M , whose atmosphere is interpolated in the grid of models 
by Allard and Haushildt (1997). Varying the convection treat- 
ment in the interior from a m = l/H p — 1 to 1.9, increases the 
temperature at the bottom of the envelope by ~ 3%. As we see, 
the density is instead increased by ~5Q%. The difference in the 
structure is sufficient to change the lithium depletion by a fac- 
tor ~10. As shown in Fig. [2 the lithium burning phase is com- 
pleted when the age is still as low as 10 7 yr. These models do not 
include mechanisms for depletion other than plain convective 
mixing. These latter, however, act on a much longer timescale 
during the MS phase, so that the remaining lithium abundance 
should be preserved at least until the age of the Pleiades (see, 
e.g. Pinsonneault et al. 1990). Notice that lithium burning be- 
gins approximately when the radiative core develops (in Fig.^ 
where the central density curve detaches from the density at the 
bottom of convection). This occurs at the knee in the HR dia- 
gram, where the track shifts from its convective behavior (the 
Hayashi track) to the so called Henyey -radiative- line of ap- 
proach to the MS (e.g. Stahler 1988). Thus the lithium burning 
begins at a T e ff around the minimum luminosity of the track, 
which we tentatively use as a possible interesting parameter. 
This r e ff_i; m gives us an indication of the PMS track location 
in the HR diagram (being smaller for a smaller efficiency of 
convection) and we expect that it might be correlated with the 
PMS lithium burning. 
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Fig. 1. Along two tracks for the solar mass, differing in a m , we 
show the lithium depletion (upper curves, scale on the right) 
and the central density and the density at the bottom of the 
convective envelope (curves which detach at log t ~ 6, scale on 
the left). 



Table 2. Lithium depletion in NextGen based models with a = 
1 and a — 1 .9 in the interior. 



M/M 






log r cff _ lim 


log N(Li) 


1.1 


1.9 


3 


3.663 


2.93 


1.1 


1.9 


10 


3.653 


2.93 


1.1 


1.9 


100 


3.650 


2.93 


1.1 


1 


3 


3.619 


3.27 


1.1 


1 


10 


3.615 


3.27 


1.1 


1 


100 


3.613 


3.27 


1 


1.9 


3 


3.649 


1.96 


1 


1.9 


10 


3.637 


2.52 


1 


1.9 


100 


3.633 


2.62 


1 


1.0 


3 


3.604 


3.111 


1 


1.0 


10 


3.609 


3.249 


1 


1.0 


100 


3.607 


3.236 


0.9 


1.9 


3 


3.631 


0.845 


0.9 


1.9 


10 


3.619 


1.920 


0.9 


1.9 


100 


3.615 


2.118 


0.9 


1.0 


3 


3.597 


2.880 


0.9 


1.0 


10 


3.593 


2.980 


0.9 


1.0 


100 


3.595 


2.930 


0.8 


1.9 


3 


3.604 


0.85 


0.8 


1.9 


10 


3.609 


1.93 


0.8 


1.9 


100 


3.607 


2.12 


0.8 


1.0 


3 


3.582 


2.345 


0.7 


1.9 


3 


3.604 


-3.36 


0.7 


1.9 


100 


3.609 


-1.49 



Table[2shows the remaining lithium abundance in the solar 
tracks after the PMS depletion, in several models from paper 
I, all of which fulfilling the requirement that the solar radius 
is achieved at the solar age. The models are either computed 
by using the new Vienna ATLAS9 (Kurucz 1993, 1995, 1998) 
model atmospheres (Heiter et al. 2002) or the NextGen models 
by Allard and Hauschildt (1997, AH97). Among the ATLAS9 
models, some are computed by the Full Spectrum Turbulence 
(FST) convection description by Canuto et al. (1996) both in 
the atmosphere and in the interior. We see that depletion is very 
severe in the FST model. For the MLT models computed with 
the Vienna ATLAS9 model atmospheres (Heiter et al. 2002) 
which have a- dtm = 0.5, the solar fit is achieved by assuming 
a larger ofo in the interior computation. The larger the r p h at 
which the atmosphere ends, the larger is the c?i n required for 
the solar fit, and the smaller the lithium depletion. In fact, the 
larger is r p h, the larger is the fraction of the over-adiabatic re- 
gion integrated using a low efficiency of convection, so that the 
r e ff of the model is significantly reduced, as well as the temper- 
ature at the bottom of the convective envelope at the epoch of 
lithium burning. Nevertheless, the minimum depletion achieved 
in these models is ~ 0.6 dex, still too much to be compati- 
ble with the open cluster LDPs (see later). Thus we confirm 
that these models using non-grey boundary conditions have a 
strong lithium depletion, as it is obtained in the grey MLT mod- 
els (e.g. Schlattl and Weiss 1999, Piau and Turck-Chieze 2002) 
as well as in the grey FST models (D'Antona and Mazzitelli 
1994, 1997). 

Table 2 shows the lithium depletion for masses from 1 . 1 M 
to 0.7 M in two sets of MLT models computed using Allard 



and Hauschildt (1997) NextGen model atmospheres as bound- 
ary condition. All these models consider a atm = 1, the value 
fixed in the grid of model atmospheres, but one set is computed 
assuming a m = 1 .9 in the interior, so that the solar track fits the 
solar radius, and in the other set we adopt the same value than 
in the atmosphere, a m = 1 , even if the corresponding evolution- 
ary track does not fit the Sun. These latter models, whichever 
is the matching point r p h, have very similar lithium depletion 2 . 

In Table |2] we also list the parameter r e ff_ii m , and Fig. |2] 
shows the lithium depletion versus r e ff_ii m : for large enough 
masses (M > 1 . 1 M ) the depletion is small for any convection 
efficiency, but the location of the track increases with the 
convection efficiency (Z" e ff_H m increases). For the lowest masses 
(0.7 M ) the r e jf location does not vary significantly, but the 
lithium depletion varies by several orders of magnitude from 
the models with low to those with high convection efficiency. 
The masses O.8-1M show a wider range of variability, and a 
very significant degree of correlation between r e ff_ii m and the 
lithium depletion. 

Fig.|2]also shows the points corresponding to the ATLAS9 
MLT models fitting the Sun: these have very similar lithium 
depletion, for the solar track, to that of the AH97 models also 

2 If the interior and atmospheric physics were exactly the same, we 
should obtain the same model structure, at least for the r p h = 10 and 
100 models (at T ph = 3 the diffusion approximation is not yet strictly 
valid and this may give some differences). Actually, also numerical 
differences can be responsible for some variation in the results: in fact 
the atmospheric structure in the AH97 model atmospheres includes 
only very few mesh points between t = 10 and 100, while the interior 
integration is much more detailed. 
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fitting the Sun. For the other masses, the results differ, but they 
follow the same pattern along the proper relation lithium - 
^eff-iim- We did not consider ATLAS9 models not fitting the 
Sun, but, if we had assumed such a low efficiency of convec- 
tion also in the interior, we would have obtained smaller lithium 
depletions, as shown by the AH97 models with a = 1 both in 
the atmosphere and in the interior. We immediately recognize 
that the lithium depletion in young open clusters requires the 
minimum efficiency of convection among those we have con- 
sidered: only the models with AH97 model atmospheres and 
<2atm = ffin = 1 can fit the open cluster lithium data, although 
the tracks do not fit the solar location. This is shown in Fig. 
[5] where we compare the lithium depletion predicted by the 
models in Table 2 with the Pleiades data by Soderblom et al. 
(1993) and Garcia Lopez et al. (1994). We have normalized 
the lithium abundances in the models to an initial abundance 
log N(Li)=3.1. We see that only upper squares, corresponding 
to the models with a dtm — a m = 1, are compatible with the 
data. In addition, all the models of O.7M have depletion larger 
than allowed by the data, suggesting that some additional phe- 
nomenon is taking place at lower T e g. In the figure we also 
plot (full line) the depletion from the models by Ventura et al. 
(1998b) computed including the thermal effect of a dynamo in- 
duced magnetic field on the convective temperature gradients. 
In this latter case, the magnetic field acts to reduce the convec- 
tive efficiency, allowing the fit with the cluster LDP. We con- 
clude that PMS depletion is compatible only with models in 
which PMS convection is much less efficient than MS convec- 
tion. 

3. Constraints on the PMS location of tracks 

We have mentioned that the (scarce) data on PMS masses seem 
to indicate that the tracks most compatible, in terms of T e s and 
luminosity, with the observations are those in which convection 
is less efficient (Covino et al. 2001, Steffen et al. 2001, Simon 
et al. 2000). A valid addition to the analysis can be the consid- 
eration of the Lithium content. An example of how the analy- 
sis could be done is shown by means of the only relevant ob- 
servation available, concerning the secondary of the eclipsing 
double-lined spectroscopic binary, containing both components 
in PMS, RXJ 0529.4+0041 (Covino et al. 2000). This system 
provides powerful constraints: it is a very important indicator 
of the fundamental stellar parameters, as masses and radii of 
the components are simultaneously derived. Fig. 3 shows the 
location of the secondary of RXJ 0529.4+0041 in the HR dia- 
gram, compared with the tracks having convection parameters 
fitting the Sun (upper panel) or not fitting the Sun (lower panel). 
The mass of this star is given as 0. 925+0. 005M o in Covino et 
al. (2001). We see that both sets of tracks are compatible with 
this mass within the errors. The tracks of 1 and O.9M in the 
sets not fitting the Sun very nicely bracket the star location and 
seem most adequate. Notice that in the following we rely on 
the assumption that both the Covino et al. 2001 data and the re- 
sults concerning the HR diagram location of PMS objects can 
be taken at face value. Actually, the errors on the T e ff determi- 
nation might result to be larger than given by the authors, and 
for this reason we judge that a more quantitative analysis is still 
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Fig. 2. Lithium versus T'efr-Um relation for tracks from 1.1 to 
O.7M , as labeled in the figure. The full dots are the mod- 
els computed with AH97 model atmospheres, while the open 
squares employ the ATLAS9 MLT Vienna model atmospheres 
by Heiter et al. (2002). The O.8M models of a- dtm = a m =\ co- 
incide at the top of the sequence of O.8M models, while they 
differ slightly for the other masses, due to possible differences 
in the input physics of the model atmosphere and interior. 

premature. The lithium abundance determination by Covino et 
al. (2000), log N(Li)=2.4±0.5 is shown in Fig. fj As the star 
may still be at a phase in which PMS depletion is going on, this 
determination is an upper limit to the future location of the star 
in the [lithium vs. MS-r e ff] plane. According to the authors, 
actually the abundance determination should be regarded as a 
lower limit, as it was done by using the Pavlenko and Magazzii 
(1996) model atmospheres for a gravity log g=4.5, more ap- 
propriate for main sequence gravities. Reduction of the gravity 
would thus provide a larger abundance. In spite of its large er- 
ror, this value, considered as a lower limit, is therefore well 
compatible with the quite low depletion of the a m = 1 models 
and with the average lithium depletion expected from the ob- 
servations in young clusters. Further studies of this system, and 
in particular a better determination of its lithium abundance, 
are important to help us in discriminating how inefficient con- 
vection must be to be compatible with these data. 

4. Conclusions 

Although the parametrizations used for convection are highly 
unsatisfactory, we have shown that there are two different in- 
dependent indications which imply that convection in the PMS 
phase must be highly inefficient: 

1 . the location in the HR diagram of the few objects for which 
the mass (or an upper value of the mass) is known, accord- 
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Fig. 3. The Pleiades data by Soderblomet al. (1993) and Garcia 
Lopez et al. (1994) (open squares) are compared with the deple- 
tion predicted by the models in Table 2 (full big squares). The 
models are placed at the r e jf they would take in an empirical 
MS, at the Pleiades age. Only the upper squares, correspond- 
ing to the models with a atm = 0[„ = 1, are compatible with 
the data. The full line shows the depletion from the models by 
Ventura et al. (1998b) computed including the thermal effect 
of a magnetic field on the convective temperature gradients. 
The large open square with the error bar represents the lithium 
abundance of the secondary component of RXJ 0529.4+0041 
(log N(Li)=2.4+0.5, Covino et al. 2001). The T eB at which the 
point is located (5500K) is assumed to be the main sequence 
T e fi of a star of mass O.925M . 



ing to several comparisons appearing in the most recent as- 
tronomical literature (although we have still to remind that 
further work is necessary to understand carefully the T e g 
scale of PMS objects). 
2. the PMS depletion of lithium, according to the well known 
observations of LDPs in young open clusters, plus the de- 
termination of the lithium abundance in the secondary com- 
ponent of the system RXJ 0529.4+0041, whose mass is 
known with high precision. 

We have also shown that the PMS location of the tracks in the 
HR diagram and the corresponding lithium depletion during 
the PMS stage are correlated, for those masses in which lithium 
depletion is significant but not complete. Any effort to increase 
our knowledge of the PMS masses and/or of the lithium content 
in young stars will add further weight to our main conclusion, 
which is that we can not parametrize the efficiency of convec- 
tion in the same way in the MS and in the PMS. A potential 
implication of this result, preliminary explored in some previ- 
ous work (Ventura et al. 1998b, D'Antona et al. 2000) is that 
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Fig. 4. The top figure shows the tracks computed with AH97 
atmospheres and ff; n = 1.9, while the bottom figure shows 
the tracks with a m = 1.0. For each mass (1.1, 1.0, 0.9, 0.8 
and 0.7M Q , from top to bottom) the three curves represent the 
tracks obtained for different r p h. The solar location is shown: it 
is compatible with the solar model having = 1 .9, but ~400K 
hotter than the solar model with ar; n = 1 .0. The location of the 
secondary component of the binary RXJ 0529.4+0041 is also 
shown. 

there might be an additional parameter playing a role in this 
game, connected with the fast rotation of the PMS stars. 
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